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Missing bits of the solar jigsaw puzzle: small-scale, kinetic 
effects in coronal studies 
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ABSTRACT 

The solar corona, anomalously hot outer atmosphere of the Sun, is traditionally de- 
scribed by magnetohydrodynamic, fluid-like approach. Here we review some recent 
developments when, instead, a full kinetic description is used. It is shown that some 
of the main unsolved problems of solar physics, such as coronal heating and solar flare 
particle acceleration can be viewed in a new light when the small-scale, kinetic plasma 
description methods are used. 



1 INTRODUCTION 



Solar corona, a tenuous and very hot outer part of the solar 
atmosphere which sits on top of other layers (photosphere, 
chromosphere, and the transition region) can be seen either 
in white light observations during the solar eclipses (because 
of 10 6 times more flux of visible photospheric photons is 
otherwise outshining it) or in extreme ultra-violet and X- 
ray observations from space (because of atmospheric absorp- 
tion o f the short wavelength radiation). ICampbell fc Moord 
describe the first successful observations of the green 
coronal line A = 5303A by W. Harkness and independently 
by C.A. Young during the total solar eclipse of August 7, 
1869. This "new" line was attributed to a new chemical 
element "coronium". It was easy to make such a mistake 
because just a year before. 1868, during the solar eclipse 
iFrankland fc Lockverl (|l868l ) discovered a prominent yellow 
line A = 5880A, which could not be ascribed to any known 
chemical element at the time. This marked a discovery of 
Helium after Greek word 'Helios' meaning 'Sun'. On Earth 
helium was only found about 10 years later by W. Ramsay. 
The blunder with coroniu m went on f or qu ite some time. 
Moreover, as described by iBirkelaridl (Il908l ). the new gas 
geo-coronium extends from 200 km above Earth into the 
space, and that the steady auroral arcs are due to electric 
luminescence in this gas. The truth had to wait until the be- 
ginning of 1940s. In hi s George Darwin lecture of the Royal 
Astronomical Society (|Edlen. B.lll945l ) entitled "The iden- 
tification of the coronal lines" Edlen describes the process 
of realising that the mysterious coronal lines, including the 
green coronal line A = 5303A , come from so-called forbidden 
lines of highly ionised metals. For example, the most intense 
green coronal line comes from FeXIV, iron atom with 13 
electrons stripped off. Edlen rightfully acknowledges that the 
correct identification of the coronal lines was triggered by a 
letter of 13 February 1937 sent to him by W. Grotrian. These 
forbidden lines are emitted from energetically unfavourable 
metastable atomic energy levels. Under normal laboratory 
conditions electrons do not readily make transitions from 
the metastable states to lower energy states (whilst emitting 



photons) due collisional de-excitation. In the extremely rar- 
efied conditions of the solar corona such transitions become 
possible. The highly ionised metal ions can only be created 
in extremely hot conditions T > few 10 6 K. Thus, at the 
time it was rather difficult to accept that solar corona is so 
hot. Especially because this seemed in an apparent viola- 
tion of the second law of thermodynamics (in Rudolf Clau- 
sius's formulation) that the heat cannot flow from a cold 
body (photosphere at « 6000 K) to a hotter body (corona 
at few 10 6 K). In this manner, one of the major unsolved 
problem in solar physics has emerged: which physical mech- 
anisms are responsible for heating of solar corona? Some 
significant progress has been made, partcitlucalry, with the 
advent of space era i n solving the coron al heating problem 
(see e.g. chapter 9 in I Aschwandenl (|2006l )). As it happens in 
science, in general, the progress has led to yet more, pre- 
viously unknown challenges such as finding the answers to 
the questions: (i) Why solar wind speed in the polar regions 
is twice higher during the solar minimum than plausible so- 
lar wind models can explain? (ii) Why the solar flares, the 
violent events in the corona when magnetic energy is con- 
verted into other forms, occur on timescales much shorter 
than the plausible resistive timescales predict? (iii) Why the 
most of the solar flare energy is taken away by the acceler- 
ated, super-thermal particles? and so on. It is important to 
realise that reproducing or mimicking the reality in (mostly) 
numerical models does not necessarily mean that we have an 
understanding of physical mechanisms in action. Here are a 
few examples e.g. in simple ID coronal loop models a steady 
state of few 10 6 K plasma can be maintained by balancing an 
ad hoc heat input with realistic losses such as heat conduc- 
tion, radiatio n, etc, or even time-va rying flare dynamics can 
be simulated (Tsiklauri et al.ll2004l). More realist ic 3D mod- 
els (|Gudiksen fc Nordlundll2005l : lGudiksenll2009t ) mimic the 
solar coronal observations and yet use unrealistically high 
magnetic resistivity. No wonder if one puts sufficient heat- 
ing or artificially adjusts dissipation coefficients, a few 10 6 
K plasma can be obtained that "looks like" the solar corona. 
But did we learn what generates this heat or makes the dis- 
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sipation coefficients anomalously high? Albeit, the answer 
is no. Notwithstanding, something useful can be learnt e.g. 
where spatially heat needs to be deposited to reproduce the 
observations. Similarly, consider the sophisticated solar wind 
models which solve numerically flui d-like equation s for elec- 
trons, protons and alpha particles l)Li et al.ll200fJ ). Indeed, 
such models can reproduce twice as high solar wind speeds 
in the polar regions than near the solar equatorial plane, as 
observed during solar minima. But, this is obtained based 
on an ad hoc energy flux injected into the both ion species. 
Again, it comes as no surprise that ad hoc energy flux (es- 
sentially an additional momentum added) to the ion species 
produces observed fast solar wind speeds. Again, no answer 
is provided to the main question: what provides this addi- 
tional energy flux (momentum)? Is it due to absorption of 
waves? or due to some microscopic reconnection events at 
the base of the corona? As to the question of why the most 
of the solar flare energy is taken away by the accelerated 
particles, the situ ation is analogous. Kinetic- scale Particle- 
In-Cell modelling (jTsiklauri fc Haruki|[2008ah can m princi- 
ple reproduce the solar flare observations in that within one 
Alfven time, somewhat less than half (^ 40%) of the initial 
total (roughly magnetic) energy is converted into the kinetic 
energy of electrons, and somewhat more than half (^ 60%) 
into kinetic energy of ions (similar to solar flare observa- 
tions). Also, a sizable fraction (up to 20%) of the magnetic 
energy can be released/conversted into other forms such as 
kinetic energy of super-thermal particles and to much small 
extent waves. However, such simulations are performed usu- 
ally in a domain sizes of several hundreds of electron De- 
bye length Ad = vt e /cj pe = 2.2 x 10~ 3 m, which is about 
10 5 — 10 6 times smaller than the s olar flare particle acceler- 
ation site. Here Vte ~ \l kBT/m e is electron thermal speed 

taken at T = 10 6 K and UJpe = yj ne 2 /(m e eo) is electron 
plasma frequency taken at n = 10 15 m -3 , commensurate 
to solar corona. The question is whether realistic behaviour 
mimicked in small simulation domains will also hold if they 
are up-scaled to realistic sizes? The situation is analogous 
to the stability of plasma in thermonuclear fusion devices 
such as Tokamaks: small-scale ("laboratory" -scale) versions 
behave quite differently from the plasma stability point of 
view than the large, "factory" -scale ones. 



2 MAGNETOHYDRODYNAMIC (MHD) VS 
KINETIC PLASMA DESCRIPTION 

The above discussion about the spatial scales, L, - the char- 
acteristic length scales of a physical system over which prop- 
erties of physical quantities (e.g. magnetic field of an active 
region in the solar corona) change - brings us to the follow- 
ing dichotomy or a dilemma. One one hand, solar corona is 
well described by the MHD approach. This is because if we 
take L as the hydrostatic scale-height of the corona (height 
over which pressure drops e-times) ~ 50 x 10 6 m, this is 
much larger than all relevant kinetic scales: (i) cjio V i ~ 7.2 
m, ion (proton) inertial length (also called ion skin depth), 
(ii) c/uipe « 0.2 m electron inertial length (also called elec- 
tron skin depth), (hi) rL, e = vte/^ce = 2.2 x 10 _3 m electron 
Larmor radius. (ui ce — eB/m e is the electron cyclotron fre- 
quency, which quantifies rotation of electrons around the 
magnetic field on their helical path.) On the other hand, 
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Figure 1. A sketch of typical power spectrum (energy per 
wavenumber fc, (i.e. E/k), ve rsus k) of magnetic fluctuat ions in 
the solar wind. Adopted from [Goldstein fc Roberts! lll995T ). 



all interesting physical processes such as wave dissipation, 
nano-scale fast reconnection (coronal plasma heating) and 
particle acceleration (during flares) occur at small, kinetic 
scales. The dichotomy, depicted in Fig. 1, is in that the bulk 
plasma dynamics in the solar corona (where most of the 
energy is stored) is observable remotely from Earth (typ- 
ical smallest coronal structures resolved e.g. with TRACE 
satellite is 0.5" i.e. « 0.36 x 10 6 m on the Sun). But the 
above interesting physical processes all operate at scales be- 
low 10 m. Thus, unless an in-situ probes are sent to the 
solar corona the kinetic-scales processes will never be ob- 
served. It seems unlikely humankind will ever come up with 
materials (for the probes) which will withstand the heat of 
the coronal environment. However, ESA has plans for Proba 
mission ( |http: / / www.esa.int /esaMI /Proba/ 1 , in which the 
telescope of Proba-3's solar coronagraph will be mounted 
on one of the spacecraft, while the other craft is manoeu- 
vred to accurately occult (mask) the main disc of the Sun. 
With the proposed Proba-3 arrangement, it is anticipated 
that accurate measurements will be possible from 1.05 to 
3.2 solar radii. Note that the discussion about spatial scales 
L can always be cast into temporal scales, T or frequencies 
f — 1 /T by using characteristic speed of processes V c h ar via 
T = L/Vchar- For example V c har can be the Alfven speed, 
Va = Bo I \/ fionrrii = Bo /^/fiop, if MHD-scale processes are 
considered. This would equally be applicable to waves and 
reconnection. In the latter case reconnection outflow speeds 
would be implied, while in the case of former - phase speed 
of the wave. 

Let us briefly discuss conditions of when the above men- 
tioned kinetic scales become important. Strictly speaking 
MHD equations are scale-free, in a sense that there is no 
physical spatial scale appearing in the equations (except L 
which is an arbitrary size of the system, which simply needs 
to stay macroscopic so that continuum mechanical (fluid- 
like) description remains valid). The devil is in the detail, 
however. In order to obtain closed set of equations one needs 
to express the electric field by means of other physical quan- 
tities. In MHD, the induction equation for the magnetic field 
is obtained by plugging in E into the Faraday's law (from 
the Maxwell's equations) dB jdt — V x E. Intrinsic physical 
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scales start appearing when instead of usua l MHD-version 
generalised Ohm's law (|Drake fc Shay||2007l ). p. 89, 

E = —v x B + nj (1) 

one starts to use a more general form: 

E — —v X B + T]j 

ne ne A \at \ ne J J ne 

Eq.(2) essentially can be obtained from the electron equation 
of motion if electron speed v e is replaced by v e — Vi — j /ne. 
The latter is mathematical statement of: when electron and 
ion speeds are different, this causes charge separation elec- 
tric field and hence current j = en(vi — v e ). In MHD, how- 
ever, ion and electron sepa rate dynamics is ig nored and 
Vi — v e = v. It can be shown (IDrake fc ShavH2007r ) that each 
of the non-MHD terms (second line of Eq.(2)) have intrinsic 
spatial scales associated with them: (i) the term containing 
j x B (called Hall term) becomes important on scales com- 
parable to c/ujpi, ion skin depth; (ii) the term containing 
(m e /ne 2 )dj /dt (called electron inertia term) becomes im- 
portant on, c/uj p e, electron skin depth scale; and finally (iii) 
the term containing the divergence of the pressure tensor, 
V • P e , (called pressure tensor term) becomes important on 
scales comparable to TL,e, electron Larmor radius. 



3 SMALL SCALE EFFECTS IN 
RECONNECTION 

In order to better understand ordering of different spatial 
scales let us consider a simple problem of steady recon- 
nection when oppositely directed magnetic field lines are 
brought into a diffusion region by an inflow of plasma, then 
they change connectivity and plasma outflow carries the field 
lines away. The change of connectivity occurs in the dif- 
fusion region which has width 5 and length A. In MHD, 
in the simplest possible formulation (Sweet-Parker model): 
(i) Bernoulli equation determines the outflow speed being 
the Alfven speed, V ou t = Va\ (ii) continuity equation pre- 
scribes in the reconnection inflow speed Vm = (<5/A)Va; 
while (iii) the generalised Ohm's law in which the advection 
term v x B is balanced by the resistive term nj sets the re- 
connection rate M sp = V in /V ut = Vi n /V A = S ,_1/2 << 1, 
where S » 1 is the magnetic Lundquist number. Thus, 
in MHD Sweet-Parker model (in which A is fixed at L) is 
producing reconnection rates much smaller what is observed 
in e.g. flares in solar corona. For the parameters commensu- 
rate to solar corona (n = 1.0 x 10 15 m" 3 , T = 1.0 X 10 6 K, 
L = 10 5 m, B = 0.01T (100 Gauss and hence V A = 6.9 x 10 B 
m s _1 ), S — 3.7 x 10 11 ) the classical Sweet-Parker rate is 
M 3 p = 1.6 x 10 -6 . The reconnection rate is also interpreted 
as the ratio of Alfven time (ta = L/Va ~ 0.0145 s) and resis- 
tive (or reconnection) times. This means in the Sweet-Parker 
model resistive (or reconnection) time is 0.0145/S _1/ ^ 2 s 
= 0.1 days. On contrary, in the observations flares last for 
up to tens of minutes. Petschek model alleviates this prob- 
lem by making the diffusion region length, A « L. This 
means that the plasma inflow speed Vi n ~ (5/A)Va a siz- 
able fraction of the Alfven speed making reconnection rates 
faster and hence flare times shorter commensurate to the 




Figure 2. A sketch of kinetic-scale magnetic reconnection. Blue 
and red curves in the quadrants show an opposite sign magnetic 
flux. The outer, blue box shows a region where ions become de- 
magnetised, whilst electrons are still magnetised. From this scale 
(c/ujpi) and below the magnetic field is carried by electrons. Black 
box with width c/LU pe shows electron diffusion region where mag- 
netic connectivity change take place. The green box with width 
shows region where electron meandering motion makes pres- 
sure to become a tensor with off-diagonal (non-gyrotropic) terms 
playing a major role. Blue and black dashed lines with arrow 
heads show ion and electron flow, respectively. Note how ion flow 
deflects at the scale c/u) p i, while electron flow does so at c/uj pe . 
The green wiggle in the centre show an electron orbit. Its radius is 
small and follows a magnetic field line away from the diffusion re- 
gion. The radius becomes large and motion chaotic (meandering) 
as the magnetic field drops to zero close to the magnetic null. 



observations. However, Petschek model is a phenomenolog- 
ical one (not derived from the first principles) and for it 
to work (as demonstrated by numerous numerical simula- 
tions) it requires plasma resistivity to be non-uniform and 
be enhanced on the diffusion region (so-called anomalous 
resistivity). Whilst, the latter being quite plausible, given 
that turbulent transport is known to enhance the dissipa- 
tion/resistive coefficients dramatically, we have not sent in- 
situ probes to the solar corona yet, and hence we do not 
know for certain whether it is turbulent or not! 

It should be noted that in the laboratory plasma Mag- 
netic reconnection experiment (MRX) measurements show 
indication of the reconnection rates or equivalently anoma- 
lous resistivity values 10s and even 100 time s large r that 
Sweet- Parker model predicts (lYamada et al.l [2006;). The 
MRX experiment plasma is not in the turbulent state. What 
is observed, however, is that the reconnection rates or equiv- 
alently anomalous resistivity values are attained when the 
transition from collisional (resistive MHD) to collisionless 
(kinetic) regime takes place. 

So, is there an alternative way (i) to have fast reconnec- 
tion i.e. get solar flare times right? (ii) also to demonstrate 
that the large proportion of flare energy goes into the accel- 
erated particles? and (iii) not to make assumptions about 
the turbulent state of the corona (hence postulating the 
anomalous resistivity)? It turns out that plasma kinetic ap- 
proach can provide answer "cautious yes" to all of these fun- 
damental questions. Here is how: For the solar coronal pa- 
rameters: a temperature of 1.0 x 10 6 K, Coulomb logarithm 
of 18.1, the Lundquist number (using Spitzer resistivity) is 
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3.7 x 10 11 . Here L — 10 5 m was used. One of the reasons 
for going beyond resistive MHD is comparing typical width 
of a Sweet-Parker current sheet S sp = S~ 1 ^ 2 L = 0.16 m to 
the ion skin depth. Typical scale associated with the Hall 
term in the generalised Ohm's law at which deviation from 
electron-ion coupled dynamics is observed is, c/u> P i = 7.2 m. 
Here particle number density of n = 1.0 x 10 15 m -3 has been 
used. Hence, the fact that c/(u} p i8 sp ) = 44 2> 1 justifies go- 
ing beyond single fluid resisti ve MHD approx i matio n (a sim- 
ilar conclusion is reached by lYamada et all l|2006T ) in their 
Fig. 12). In other words, ion skin depth scale is reached first 
before resistive MHD Sweet-Parker current sheet. Let us put 
this fact (that c/(u> p iS ap ) 1) in the visual context using 
Fig. 2: On a global (bulk MHD) scale L electrons and ions in- 
flow towards the diffusion region (where magnetic field con- 
nectivity changes) in a coordinated way as if they are glued 
to each other. When the flow approaches ion skin depth scale 
c/ui P i ions become demagnetised (magnetic field is no longer 
is frozen into ion fluid) and magnetic field is carried forward 
by electrons. At c/ui p i scale the Hall term in the generalised 
Ohm's law dominates. This picture is broadly corroborated 
by kinetic, particle-in-ce ll simulations of x-point collapse 
jTsiklauri fc HarukillioOTl ) . In panel (a) of Fig.3 we see that 
for electrons the diffusion region width is indeed c/uj pe , while 
for ions, panel (b) of Fig.3, the width of the region where 
ions start to deflect from the current sheet (region where ions 
become demagnetised) is about c/uo P i = 10c/cu pe (note that 
here the ion to electron mass ratio of mi/m e = 100 has been 
used). This continues until electron skin depth scale c/uj pe 
is reached where electron inertia term takes over. Further 
down the spatial scales, when electron Larmor radius scale 
rL, e is reached electron mean dering motion (dashed green 
wiggled curve in Fig.2) sets in (|Horiuchi fc Satoll 19971 ). This 
makes the non-gyrotropic (off-diagonal) components of the 
electron pressure tensor to dominate other terms in the gen- 
eralised Ohm's law, and hence makes them responsible for 
breaking the frozen-in condition. The concept of pressure 
tensor is somewhat counter-intuitive. This can be compre- 
hended as follows: away from the electron diffusion region 
pressure is scalar (isotropic) and electrons flow along (are 
frozen into) the magnetic field lines. Inside the the diffu- 
sion region ^ rL, e electrons become demagnetised and ex- 
ecute meandering (rather chaotic) motions. These are such 
that pressure becomes different as one travels in different 
directions (i.e. pressure becomes a tensor as opposed to a 
scalar). Previous results on collisionless reconnection both 
in tearing unstable Harris current sheet jK uznctsov a et al.l 
1 19981 ; iHesse et al.lll999l ; iBirn et al.ll200ll;|Pritchetdl200ll) and 
stressed X-point collapse ( Tsiklauri fc Harukil 120071 12008a) 
have shown that magnetic field is frozen into electron fluid 
and the term in the generalised Ohm's law that is responsi- 
ble for breaking the frozen-in condition is electron pressure 
tensor off-diagonal (non-gyrotropic) component gradients. 
Thus, there is a case for inclusion of the electron pressure 
tensor non-gyrotropic components in a model of co llisionless 
recon nection. Analytically, this was achieved by (|TsikIauril 
2008). In this simple model instead of balancing the ad- 
vection term in the generalised Ohm's law with the resistive 
term (as discussed above), non-gyrotropic components of the 
electron pressure tensor gradient were used. This produced 
simple analytical results which explain well (i) the recon- 
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Figure 3. Electron (a) and ion (b) flow patterns at the peak 
of time-transient re c onnec tion of x-point collapse. Adapted from 
iTsiklauri fc Harukil J2007h . who used particle-in-cell, kinetic ap- 
proach for the numerical simulation. The arrow size represents in- 
stantaneous, cell-averaged velocities normalised by speed of light. 
Note that for electrons the diffusion region width is about c/uj pe , 
while ion flow deflects from the diffusion region at a width c/u> p i 
(here c/u> p i = 10c/u> pe ), as sketched in Fig. 2. 



nection rate and (ii) width of the electron diffusion region 

measurements in MRX experiment. 

There is a growing amount o f work ( Dra ke fc Shavl2007l ; 
ITsiklauri fc Harukil I20071 l2008aT ) that suggests that in the 
collisionless regime, on the scales less than c/ui p i magnetic 
field is frozen into the electron fluid rather than bulk of 
plasma. One can write in general v e = Vi — j /{en). This 
relation clearly shows that in collisional regime (when the 
number density n is large), the difference between electron 
and ions speeds diminishes v e — Vi = v. However, as one 
enters collisionless regime (when the number density n is 
small) the deviation between el e ctron a nd ion speeds starts 
to show. In ITsiklauri fc Harukil (|2008al ) we proposed a pos- 
sible explanation why the reconnection is fast when the Hall 
term is included. Inclusion of the latter means that in the re- 
connection inflow magnetic field is frozen into electron fluid . 
As it was previously shown in ITsiklauri fc Harukil (|2007l ) 
(see their Figs. (7) and (11)) speed of electrons, during the 
reconnection peak time, is at least 4-5 times greater than 
that of ions. This means that electrons can bring in / take 
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out the magnetic field attached to them into / away from the 
diffusion region much faster than in the case of single fluid 
MHD which does not distinguish between electron-ion dy- 
namics. Thus, it is clear that inclusion magnetic field trans- 
port by electrons is crucial in resolving of above mentioned 
problems MHD description of the solar corona faces. 

To summarise, some of the time-dependent collision- 
less reconnecti on models such as tear i ng unstable Harris 
current sheet (Kuznctsovaetal ] Il998l : iHesse et all Il999l : 
iBirn et al.ll200ll: |Pritchetdl200ll ) "and s tressed X-point col- 
lapse dTsiklauri fc HarTikifl2007l . l2008af ) where able: (i) to 
show that a sizable fraction of the magnetic energy can be 
converted/released into heat and acceler ated, super-thermal 
partic le energy (see e.g. Fig. 13 from iTsiklauri fc Harukil 
(2007)); (ii) obtain fast reconnection rates and hence mag- 
netic energy release rate commensurat e to solar flare ob- 
servat ions (see e.g. Figs. 3 and 9 from ITsiklauri fc Harukil 
Hoot})); (iii) reproduce the observational fact that the large 
proportion of flare energy goes into the accelerated particles 
in the correct partition, i.e. somewhat less than half (40%) of 
the initial total (roughly magnetic) energy is converted into 
the kinetic energy of electrons, and somewhat more than 
half (60%) into kinetic energy o f ions (see e.g. Figs. 3 and 4 
from ITsiklauri fc Harukil (|2008a! )); (iv) not to invoke anoma- 
lous resistivity. However, despite this success there are main 
challenges ahead: (i) the processes described by the collision- 
less reconnection models are small scale (< c/ui P i ~ 10 m). 
Indeed, on scale less than ion skin depth magnetic field is 
advected into the diffusion region by electrons and since elec- 
trons are lighter and move faster, the reconnection is accord- 
ingly fast. But, sufficiently far away from the electron dif- 
fusion region, electron flow must slow down to bulk plasma 
speeds when MHD description takes over (ions and electrons 
become glued to each o ther), where and how t h is happens? 
Recen t PIC simulations iDaughton et alj (|2006l ) ; IShav et al.l 
(2007) suggest that the electron diffusion region length can 
extend for much longer distances downstream than previ- 
ously thought. Why does this happen is not understood, (ii) 
How this small-scale dynamics on the kinetic scale translate 
into MHD scales? i.e. despite producing fast reconnection 
rates, these processes happen in physically small volumes. 
Will they affect MHD-scale (say a L — 10 6 m) processes? 
It will be impossible to find out until we can perform a nu- 
merical simulation 10 6 mx 10 6 mx 10 6 m using realistic coro- 
nal active region fields as an input for Particle-In-Cell or 
Vlasov simulation. The latter would seem many decades 
away even if Moore's law (that the number of transistors 
that can be placed inexpensively on an integrated circuit 
has increased exponentially, doubling approximately every 
two years) continues to hold in the future. Until such times, 
we can do lOmxlOmxlOm kinetic simulation and hope for 
the best that our models will scale correctly towards the 
larger (physically meaningful) spatial scales. 



4 SMALL SCALE EFFECTS IN WAVES 

Study of MHD waves in the solar corona has seen steady 
progress over last two decades. There are two aspects to 
the study: (i) waves carry information about the medium 
they propagate in. This is used in the field of solar coro- 
nal seismology to infer physical parameters of the corona 



see fo r reviewsfNakariak ov fc Verwichte (2005); De Moortell 
20091 )). For example traditional method of measuring mag- 
netic field in solar photosphere is based on the Zeeman ef- 
fect. However, in solar corona plasma is too dilute render- 
ing this method unusable. Based on the seminal theoretical 
work of Id wm fc Roberts! (|l983l ). transverse oscillations of 
solar coronal EUV loops were used to successfu lly measure 
the mag netic field (|Nakariakov fc Ofmanl l200lh . (ii) waves 
transport mechanical (and magnetic) energy and if dissi- 
pated may produce plasma heating. 

In the MHD approximation, in a uniform plasma pene- 
trated by a uniform magnetic field, Bo, there are three dis- 
tinct types of waves: fast and slow magnetosonic and Alfven 
waves. 

Fast magnetosonic waves are compressible MHD waves 
which have an advantage of propagating across the magnetic 
field. This implies that the heat produced by the wave dis- 
sipation can spread across and cover a large area. The fast 
magnetosonic waves that are generated in the solar interior 
cannot reach the corona because of the reflection on the pres- 
sure/density gradient (due to gravitational stratification). 
Thus, initially they were excluded from coronal heating 
models. However, in the transversely inhomogeneous plasma 
there are two possibilities to couple Alfven waves (which 
have no problem of crossing the gradient) with fast magne- 
tosonic waves : either by non- linear coupling in 2.5 MHD ap - 
proximation l|Nakariakov et a l. 1997; Tsikla uri et al.ll200ll ); 
or linear (m uch more efficient) coupling in 3D MHD ap- 
proxi mation (|Tsiklauri fc Nakariakovl 120021 : ITsiklauri et al.l 
l2003ft . 

Slow magnetosonic waves are also compressible, but 
propagate preferentially along the magnetic field. Normally, 
these waves are discounted from coronal heating candidate 
list on the grounds of no t carrying (less by four orders of 
magnitude) sufficient flux lDe Moortell (2009). However, the 
latter estimate is based on the as sumption of just a sin- 
gle ha rmonic (fixed frequency) and ITsiklauri fc N akariakov 
l|200 if ) have demonstrated that if a wide spectrum (contin- 
uous spectrum with a plausible frequency (or wavenumber) 
range) of slow magnetosonic waves exists in the corona, than 
their dissipation would provide sufficient heating. Whether 
such spectrum exists is an open question. However since in 
the nature there is always a cascade of energy from large 
scales to small scales (see e.g. Fig. 1), this wide spectrum 
conjecture seems plausible. 

Alfven wave is an incompressible, transverse wave (such 
as e.g. electromagnetic wave), with frequency, uja, much 
smaller that ion cyclotron frequency, ui C i — eB/rrii. In an 
Alfven wave the background magnetic field tension provides 
restoring force, while plasma ions provide inertia (the two 
key ingredients for any oscillatory motion). Alfven waves 
are easy to excite, e.g. a sheared plasma flow across mag- 
netic field. However, they are notoriously difficult to dis- 
sipate because of small resistivity of the solar corona, i.e. 
Alfven wave is a good vehicle for transporting energy (e.g. 
from convective motions below the photosphere into the 
corona), but this "vehicle" has virtually no breaks (hence 
cannot be stopped!). Thus, in order to deposit energy in 
the first ~ 50 x 10 6 m of the corona (the hydrostatic scale- 
height) above the transition region, some enhanced dissipa- 
tion mechanisms where proposed in the past. In the solar 
corona, basic magnetic structure that is encountered is ei- 
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Figure 4. A sketch of parallel, at an angle and perpendicular 
propagation of a wave. 



ther a closed loop (during 11-year cycle minimum more of 
these are found near solar equator and low latitudes); or 
open magnetic strictures such a plumes (during the solar 
minimum these are mostly located near the north and south 
polar regions). The common feature of these coronal struc- 
tures is the density variation across the magnetic field, i.e. 
inside the loops and plumes density is enhanced by 3 - 10 
times compared to the surrounding plasma. In the case of 
closed magnetic structures this enhanced density is believed 
to be due to material evaporation from denser layers of the 
sun such as transition region and chromosphere. In the open 
magnetic structures the density enhancement is either due 
to a solar wind and again material evaporation from below 
(which could be either due to reconnection or any other pro- 
cess which increases plasma temperature so that the balance 
no longer holds). There are two mechanisms for the above 
mentioned enhanced Alfven wave dissipation and absorption 
in the solar corona: phas e mixin g dHeyvaerts fc Priesdll983h 
and resonant absorption llonso n] [|l978h . In both mechanisms 
transverse inhomogeneity of plasma (i.e. plasma with den- 
sity variation across the uniform magnetic field) plays a cru- 
cial role (albeit for di fferent reasons). B r ief rev iew of both 
topics can be found in |Pe Moortel et al.1 i|200Sl ). 

For particle acceleration and efficient non-resonant 
wave-particle interactions (e.g. converting wave energy to 
particle kinetic energy and heat if inter-particle collisions 
are efficient), it is necessary for a wave in question to have 
magnetic field-aligned electric field. The word non-resonant 
is crucial here because, interaction of waves with particles 
when (i) the wave phase speed V p h = u)/k coincides with 
particle (thermal) speed v t h,a = \l ks T/m a leads to the 
Landau resonance; or (ii) to the cyclotron resonance, when 
the rotation frequency of the electric field vector of the cir- 
cularly polarised wave, ui coincides with the particle's cy- 
clotron frequency, uj ca = eB/m a (and when both the wave 
electric field vector and particle rotate in the same direc- 
tion). In the both cases particle experiences almost con- 
stant electric field, which facilitates efficient wave-particle 
resonant interaction. Away from the resonance V p h ^> or 
<C v t h,a in the case of Landau damping; or lo 3> or <g ui ca 
in the case of cyclotron resonance, wave electric field oscil- 
lates too fast or too slow and has no effect on plasma par- 
ticles of species a. In the MHD approximation Alfven wave 



propagates strictly along the magnetic field (with some fi- 
nite parallel to the magnetic field wave number fcii 7^ 0) 
as shown in case A, Fig. 4. In this case the perpendicular 
wave number fcj_ = 0. For the field aligned electric field to 
exist k± 7^ condition should be true. Case C, Fig. 4 also 
shows why the fast magnetosonic waves are believed to be 
efficient particle accelerators. This because they preferen- 
tially propagate across the magnetic field and hence posses 
finite k± 7^ 0. For Alfvenic waves k± 7^ (see case B, 
Fig. 6) only possible when kinetic- scale effects are included . 
In particular, it has been realised (|Hasegawa fc Chenlll976h 
that when externally applied wave frequency becomes res- 
onant with the Alfven frequency uja = fcyVk, a resonance 
occurs and the driving wave converts into kinetic Alfven 
wave (KAW) which has the finite field-aligned electric field. 
More generally, generation of the parallel (magnetic field- 
aligned) electric field occurs when the applied low-frequency 
(uj < tJci) Alfven wave has k± comparable to any of the 
kinetic spatial scales. Such Alfven waves (with E\\ 7^ 0) 
are c alled dispersive Alfven waves (DAW) l|Stasiewicz et al] 
2000). The latter split into two classes: kinetic Alfven waves 
(KAW) and inertial Alfven waves (IAW). The distinction 
between the two is drawn according to the value of plasma 
beta, which is the ratio of thermal and magnetic pressures, 
i.e. (3 — p/(B 2 /2fio). Plasma beta can also be related to 
the ratio of the sound speed and Alfven speed squared, 
i.e. f3 = {2/j)(cl/V%) w c 2 /V^f, the ratio of specific heats 
7 in e.g. adiabatic case is 5/3. Here c a = \pyp~J~p is the 
sound speed. Physically, plasma beta prescribes importance 
of either thermal effects or magnetic field. E.g in solar 
corona j3 <C 1 meaning that magnetic fields play a dom- 
inant role in coronal plasma behaviour. In the context of 
dispersive Alfven waves plasma beta prescribes which term 
in the generalised Ohm 's law supports parallel electric field 
l|Stasiewicz et a l. 2000). If /3 < m e /mi, i.e. v t e,vu < va then 
the dominant mechanism for generation of parallel electric 
field in the generalised Ohm's law is the electron inertia 
and one has inertial Alfven waves (IAW). If /3 > m e /m,i, i.e. 
vte,vu > va then the dominant mechanism for generation of 
parallel electric field in the generalised Ohm's law is electron 
pressure gradient and one has kinetic Alfven waves (KAW). 

In the solar corona, particle acceleration e.g. during so- 
lar flares is generally believed to be occurring as a result of 
magnetic reconnection. For descri ption of particle dynam- 
ics e i ther test particle approa ch l|Zharkova fc Gordovskyvl 
120051 : iDalla" fc Browning! I2OO8T ) which is ignores self consis- 
tent fields or more recently kin etic, particle-in-cell simula- 
tion (jSiverskv fc Zharkovall2009l ') has been used. The prob- 
lem with particle acceleration during solar flares is associ- 
ated with the numbers of particles involved. In the stan- 
dard flare model, reconnection event occurs higher up the 
corona and as result 50-80 % of the flare energy is con- 
verted into accelerated particles (electrons and ions). These 
rush down towards denser layers of the solar atmosphere 
guided by coronal magnetic field structures and produce X- 
ray emission via bremsstrahlung. The number of detected 
X-ray photons can be related (via certain model assump- 
tions) to the number of accelerated electrons. For a typ- 
ical flare more than 10 39 electrons are accelerated, yield- 
ing a huge current sheet (with a volume of 10 28 cm 3 ) that 
must remain stable for more than 60 seconds. This is an 
unlikely course of events because of the dissolution (thin- 
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Figure 5. The left column shows three snapshots of Alfven wave magnetic field component. The middle columns shows the generated 
parallel electric field, . The right column shows parallel velocity phase space of electrons. See text for details. 



ning) of the current sheet (1 Brown Sz Melrose! H977I ). In ef- 
fect, the above numbers imply that all of the electrons in 
the current sheet need to be evacuated. Another problem 
is with the return currents. A beam of electrons injected 
into plasma is known to be compensated by the generated 
counter beam. Thus, expecting a bunch of all 10 39 electrons 
happily travels towards the footpoints of the magnetic struc- 
ture is naive. A viable alternative for transporting flare en- 
ergy by the beam of accelerated particles, is waves. As we 
know waves can transport energy without mass transport. 
Therefore, if flare energy is delivered to footpoints by waves 
then the above mentioned difficulties (the number problem 
and return currents can be avoided). In this scenario, recon- 
nection flare event produces dispersive Alfven waves (which 
posses parallel electric fields). These subsequently travel to 
the footpoints, then accelerate particles in the vicinity of 
footpoints. The accelerated particles (mostly electrons, be- 
cause these waves are known to preferentially accelerate elec- 
trons due to their small inertia, whilst producing ion heating 
- broadening of their distribution function), in turn, produce 
the observed X-ra ys via the brems s trahlu ng. Using particle- 
in-cell simulation, Tsikl auri et alj (|2005l ) have recently ex- 
plored how circularly polarised Alfven wave with u = 0.3cj C i 
propagates in a transversely inhomogeneous plasma. This 
frequency is somewhat higher than that for MHD wave and 
yet smaller that ion cyclotron resonance. Such wave strictly 
speaking is called Ion cyclotron rather than Alfven wave. 
In the uj <C u C i limit, however, is behaves like an Alfven 



wave. The density in the middle of simulation domain was 
smoothly increased by a factor of 4 (a ramp-like function 
p{y) = 1 + 3exp [—((y— 100)/50) 6 ]) across the uniform 
magnetic field in order to model a coronal loop (accord- 
ingly temperature was varied as inverse of density so that 
total pressure balance is preserved). The widths of the trans- 
verse inhomogeneity on each side of the ramp is « c/u p i. 
Fig. 5's left column demonstrates that initially plane Alfven 
wave front becomes distorted (middle part travels slower) 
because its phase speed (Va(v) = i?o/\/ fJ>op(y)) is a func- 
tion of density, which in turn is a function of a coordinate 
y across the uniform magnetic field applied along a;-axis. 
It is also interesting to note that in the strongest density 
gradient regions y = 50 and 150 the propagating Alfven 
wave damps. In the middle column of Fig. 5 parallel elec- 
tric field E\\ = E x is shown. At time t = 0, E\\ is zero 
everywhere, but as the Alfven wave propagates and its wave 
front becomes distorted E\\ starts to grow. It typically at- 
tains values of few 10 _3 (meCu> pe /e). For solar coronal con- 
ditions this is typically exceeds the Dreicer electric field 
(which is associated with the particle acceleration runaway 
regime l|Dreicerlll959l )) 10 6 times! Hence, the generated par- 
allel electric field efficiently accelerates electrons, as can be 
seen in the right column of Fig. 5. The latter shows paral- 
lel to the magnetic field phase space (V x vs. x) where each 
dot corresponds to an electron (total of ~ 5 x 10 8 ). We 
gather that in the regions where Eu is generated, the num- 
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ber of particles with increased field- aligned velocity V x is 
increased (i.e. particles are accelerated). Note that this coin- 
cides w ith the regions where Alfven wave enhanced damping 
occurs. iTsiklauri fc Harukil (|2008bh performed a parametric 
study, exploring how the new mechanism of parallel electric 
field generation and acceleration of electrons discovered by 
ITsiklauri etafl (|2005f ) depends on problem parameters, such 
as the variation of frequency and amplitude of the applied 
(driving) Alfven wave, as well as plasma beta, affect levels 
attained by the E« and the fraction of accelerated particles 
(the latter is defined as percentage of electrons with speeds 
above the electron thermal speed in the density gradient 
regions). It was established that (i) Eu is always orders of 
magnitudes greater than Dreicer electric field and (ii) the 
faction of accelerated electrons ranges within 20-50%. The 
latter makes the discovered electron acceleration mechanism 
quite efficient and potentially capable of resolving the above 
mentioned problems of standard flare models which use 
reconnection-produced electron beams instead of waves as 
the means of delivering energy flux to the footpoints in order 
to produce the X-rays. It should be noted tha t the transverse 
inhom ogeneity is crucial in the mode l of iTsiklauri et all 
(2005j) (see also further analysis paper l|Tsiklaurill2007l )V In 
fact, when the transverse density inhomogeneity is removed, 
i.e. when Alfven wave front remains always perpendicular 
to the uniform background magnetic field, no parallel elec- 
tric field gen eration or electron acc eleration is observed (see 
Figs.7-10 in ITsiklauri et al.l l|2005h ). In terms of the above 
k± discussion, this can be explained as following: When the 
density across the magnetic field is inhomogeneous, despite 
the fact that initially k± — as the Alfven wave front starts 
to deform k±_ is efficiently generated. When on the other 
hand density is constant everywhere then Alfven wave front 
is always at right angles to the magnetic field, i.e. Alfven 
wave propagates st rictly along the field , and t hus k± — at 
all times. Recently. iFletcher fc Hud son (2008) revisited the 
problem of sola r flare electron accelerat i on wi th Alfvenic 
pulses. Further, iMcClements fc Fletcher! (|2009l ) quantified 
the full details, namely, exploring how fraction of accel- 
erated particles is affected by such crucial parameter as 
transverse length-scale l± ~ 2n/kj_ of the Alfvenic pulse. 
They found that significant fraction of electrons would be 
accelerated if l± is few meters or le ss. A clear dist i nction 
should be made between models of ITsiklauri et~al] (|2005l ) 
and IMcCleme nts fc Fletcher! (|2009l ). As explained above, in 
the case of former, k± is initially zero and it seems that 
such wave is more plausible to exist in nature. This is be- 
cause it represents simple, ubiquitous Alfvenic harmonic 
wave. Recall that even in MHD a pproximation k ± is al so 
zero at all times. It is crucial in ITsiklauri et al. | i|2005h 's 
model that plasma has density inhomogeneity mimicking 
transverse density structure of the coronal loop. It is this 
inhomogeneity that makes Alfven wave front to distort and 
create finite k± and h e nce g enerate Eu . In the case of 
IMcClements fc Fletcher! i|2009r )'s model k± / (and hence 
E\\ 7^ 0) from the start and there is no transverse density 
inhomogeneity. Also, a pulse instead of harmonic have is 
considered. If we try to compare plausibility of the two mod- 
els the following considerations come to mind. Both models 
"work" i.e. generate large enough and accelerate suffi- 
cient proportion of e lectrons if the following are true: In 
ITsiklauri et~aH (|2005l Vs model a large scale Alfvenic, har- 



monic wave initially has k± — 0, interacts with the trans- 
verse density inho mogeneity which has a scale of Unhom — 
c/ujpi ~ c/ujci. In IMcClements fc Fletcherl (|2009l Vs model 
(see their Eq.(28)) Alfvenic pulse with k± ~ l/[A(c/u) p i)] 
travels in the homogeneous plasma (where A is the pertur- 
bation amplitude). In effect, both models imply either small 
scale transverse inhomogeneity of the medium or small scale 
transverse wave number. It is unclear at this stage whether 
we have such small scale phenomena in solar corona. Per- 
haps, in situ coronal probes or remote sensing from distances 
much small than Sun-Earth distance could shed some light, 
as this is currently possible in Earth magnetospheric appli- 
cations. However, until such times the jury on this topic will 
remain out. 

In summary, wave based models offer attractive solu- 
tions to the above mentioned outstanding problems in solar 
flare particle acceleration. It is yet to be clarified however, 
how small-scale (kinetic-scale) phenomena feeds into the 
large-scale (MHD-scale) phenomena. From the fundamental 
point of view, of course kinetic-scale modelling includes all 
the essential physics. However, realistically shortcomings of 
man-made even the largest parallel supercomputers will un- 
likely offer us the performance we need. E.g. a full 3D Vlasov 
code 128 grids in the 6-dimensional phase space (3 spatial 
dimensions and 3 velocity dimensions) require 32 Terabytes 
(32,000 Gb) of RAM. Given that spatial grid size is the De- 
bye length (Ad = v t h, £ /Lo P e = 2.2 x 10 -3 m in the solar 
corona), in 3D this offers rather uninterestingly small phys- 
ical volume. Thus, there is always a restriction imposed to 
consider a lower dimensional physical systems. 



5 CHALLENGES AHEAD 

Despite significant progress in understanding of physical 
processes ongoing in solar atmosphere - it is enough to take 
a look as NASA or ESA webpages describing a number of 
previous and existing successful space missions, as well fu- 
ture mission plans - there is still a long and exciting way 
to go! As far as coronal heating problem is concerned, in 
author's opinion further progress could be made if following 
considerations are noted. 

From the observational perspective it would be rather 
important to investigate heating of active regions (ARs). 
AR is a volume of solar atmosphere above sunspots. They 
are important because during solar maximum ARs can 
be responsible for m 80% of heating of the entire corona 
|Aschwanden et alj|2007h . Top left of Fig.6 shows TRACE 
171 A image of an AR. The box in Fig.6 shows potential 
extrapolation (using Green's function method) of magnetic 
field lines based on an input from the measured normal com- 
ponent of the photospheric magnetic field that can be used 
for numerical modelling. Ideally, if we have simultaneous, 
high time cadence 2D imaging and 2D Doppler shift data at 
several wavelength, A, corresponding to emission from e.g. 
white light (continuum, photosphere, log(T) = 3.7), Hell 
(A = 304A, chromosphere, transition region, log(T) = 4.7), 
FelX (A = 17lA, quiet corona, upper transition region, 
log(T) = 5.8) FeXII, FeXXIV (A = 193A, corona and hot 
flare plasma, log(T) = 6.1,7.3), etc then one would able 
to Fourier transform the data to see e.g. how much wave 
power spectrum is lost (dissipated) from one height (a hor- 
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Figure 6. Top left panel shows TRACE 171 A image of an active 
region in solar corona. The box shows potential extrapolation (us- 
ing Green's function method) of magnetic field lines based on an 
input from the measured normal component of the photosphcric 
magnetic field that can be used for numerical modelling. 



izontal cross-section through the box in Fig. 6) to another. 
Essentially, different line forming ions probe different tem- 
peratures and therefore different heights in the box in Fig. 6. 
This would ultimately enable us to make a judgement to 
what extent waves contribute to the AR and coronal heat- 
ing? Also, measuring the abundance of the observed bi-polar 
jets, which are believed to be signatures of small-scale re- 
connection, could enable us to make a judgement on recon- 
nection^ contribution to the coronal heating. NASA's At- 
mospheric Imaging Assembly (AIA) on board of the Solar 
Dynamics Observatory (SDO), to be launched in Novem- 
ber 2009, will make such unprecedented progress. However, 
AIA's resolution of about 1 arcsec and time cadence of 10 
seconds or better will naturally limit the types of waves it 
can detect. Perhaps a better job could be done by a ground 
based project called Rapid Oscillations in Solar Atmosphere 
(ROSA), http://star.pst.qub.ac.uk/rosa/, which will have 
0.1 arcsec resolution and will take in 30-125 frames per sec- 
ond! However, ROSA being a ground based instrument, can 
only probe photosphere/chromosphere/TR lines, i.e. cannot 
probe corona, and yet a whole wealth of new data is expected 
which will hopefully lead to new discoveries. On a negative 
note, with SDO's AIA and ROSA importance of kinetic ef- 
fects still cannot be probed, because 1 arcsec on the Sun is 
0.72 x 10 6 m - orders of magnitudes larger than any kinetic 
scale discussed above. The concept of measu ring wave dissi- 
ation with height has been used recently bv lBaneriee et all 
2009). This was based on variations in EUV line widths of 
coronal plumes (viewed off the solar limb) with height. 

Also, there seems to be is a serious gap in the stud- 
ies of the coronal heating problem, as here are no stud- 
ies which consistently would monitor global (average) coro- 
nal temperature. There are studies that monitor TSI (to- 
tal solar irradiance) and its correlati on to th e (11 year-) 
solar activity cycle. Recent study by iFrohlichl l|2009h sug- 
gests that the long-term trend of TSI is most probably 
caused by a global temperature change of the Sun (pho- 
tosphere) that does not influence the UV irradiance in the 
same way as the surface magneti c fields . There is however, 
a work of lAschwanden fc Actor] (|200ll ). where they split 



the differential emission measure (DEM) into 10° sectors 
and estimate coronal heating requirement sector-by-sector. 
For every full-disk dataset with EUV and SXR coverage, it 
would be straightforward to compute a DEM of the entire 
corona, from which an emission measure-weighted temper- 
ature could be extracted as function of time. i.e. it could 
be possible to produce a graph of averaged (by all sectors 
i.e. full disk) coronal temperature as function of time, over 
the time interval of e.g. 3 x 11 = 33 years with sufficient 
time cadence (e.g. once per year or more). Such plot would 
answer a significant question: whether the heat release in 
the corona is indeed related to the magnetic field (which 
we know changes on 11 year timescale). There is certainly 
a good correl ation between the magnetic flux and SXR in- 
te nsity (e.g. (|Benevolenskava et al.ll2002j ). or see Fig. 1.13 
in lAschwanden! (120061 ) - Indeed, there seems to be a good 
case for IsxTif) cc B 2 ). However, it is questionable whether 
IsxT(t) is a good proxy for global average coronal tempera- 
ture, because SXR emission is no n thermal (produced by the 
forbidden lines) - cf. Fig. 2. 3 from I Aschwanden! (|2006t ) show- 
ing strong deviation of coronal emission from the blackbody 
spectrum. The analogy with a patient who has caught cold is 
relevant here: if a patient is sick, its temperature is routinely 
monitored by doctors. If we claim there is a coronal heating 
problem, indeed, solar physicists should monitor the global 
average temperature of their patient - the Sun's corona! 

As far as modelling challenges are concerned, it seems a 
real progress can be made if two-fluid simulation of heating 
release in AR is attempted. When ion and electron dynam- 
ics is decoupled, we know that reconnection will proceed 
fast jKuznetsova et al.lll998l: Hesse et aT] 1 19991: iBirn et al.l 
l200ll : |Pritchettll200ll : lTsiklauri fc Harukill2007l . l2008al ) even 
without invoking anomalous resistivity. What is encourag- 
ing is that we do not need to use realistic ion to electron 
mass ration of rrii/m e = 1836. As previous results have 
shown, as long as ion-electron dynamics is decoupled sim- 
ulation with mi/m e = 1 00 is as good IS esse et aT] l200ll : 
iTsiklauri fc Haruki|[2008ah . This implies that significant re- 
duction in CPU requirements can be achieved. The two- fluid 
simulation results then can be compared with the SDO's 
AIA and ROSA data to see whether inclusion of the two- 
fluid effects can provide fast enough reconnection or gener- 
ated wave power and hence adequate coronal heating. 
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